The Great Basin Naturalist 


PUBLISHED AT PROVO, UTAIE, BY 
BRIGHAM YOUNG UNIVERSITY 


ISSN 0017-3614 


VOLUME 50 30 JUNE 1990 No? 2 


MAYFLY GROWTH AND POPULATION DENSITY IN 
CONSTANT AND VARIABLE TEMPERATURE REGIMES 


Russell B. Rader’? and James V. Ward! 


ABSTRACT.—The thermal equilibrium hypothesis predicts that aquatic insect body size/fecundity and, consequently, 
population density and biomass will be maximized in geographic areas or along altitudinal gradients where the thermal 
regime is optimal with respect to growth and development. Seasonal growth analyses of three maytly species, 
combined with detailed thermal descriptions, were used to explore differences in body size and fecundity at three sites 
with similar elevations but different temperature regimes. Site 1 was located near the upper altitudinal distribution for 
each species, whereas sites 2 and 3 were located beiow a deep-release storage reservoir. The temperature pattern at 
site 1 had rapid seasonal changes, with a short summer and a long, freezing winter. Site 2 demonstrated gradual 
seasonal changes combined with winter warm and summer cool temperatures. Site 3 was intermediate with respect to 
seasonal change and winter harshness but had the highest maximum and mean annual temperatures. Mayfly develop- 
ment at site | was characterized by slow growth during the summer-autumn period, no growth during the winter, and 
a rapid increase during the spring-summer period. In contrast, growth at site 2 was continuous throughout the year, 
including the winter. Growth at site 3 was either continuous across sites or rapid during the spring-summer period, 
depending on the species. Based upon the interactions among temperature, body size, and metabolic costs, the 
thermal equilibrium hypothesis was successful at predicting body size and fecundity differences among sites. It was less 
successful at predicting variation in population density and biomass. Density-dependent and density-independent 
sources of mortality, including temperature, may interrupt the translation of higher fecundity into higher population 
density and biomass. 


Temperature, because of its influence on optimal with respect to growth, development, 
metabolism, growth, and reproductive suc- and body size. Fecundity, an essential compo- 
cess, is a dominant ecological determinant of — nent, but not the only factor defining repro- 
the geographical and altitudinal distributions ductive success, should decrease in tempera- 
of aquatic insects (e.g., Vannote and Sweeney — ture regimes warmer or cooler than optimal. 
1980, Ward and Stanford 1982, Sweeney Other factors, which may influence fecundity 
1984, Ward 1986). The thermal equilibrium — and may or may not be influenced by temper- 
hypothesis (TE) is a conceptual model of the — ature, can determine population size and dis- 
effects of temperature on aquatic insect tribution (e.g., egg-hatching success, emer- 
metabolism, growth, body size, and therefore — gence success, mating success, feeding rates, 
fecundity (Sweeney and Vannote 1978, Van- assimilation efficiency, food quality and quan- 
note and Sweeney 1980). It predicts that pop- tity, biotic interactions). The TE hypothesis, 
ulation density, distribution, and stability however, attempts to define the influence of 
(Connell and Sousa 1953) are determined by temperature on population size and distribu- 
individual reproductive success and will be — tion based only on the effect of temperature 
maximized in geographic areas or along altitu- on metabolism, growth, and therefore body 
dinal gradients where the thermal regime is — size/fecundity. 
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The objective of this study was to analyze 
the influenee of temperature on the growth 
and body size/fecundity of three maytly spe- 
cies and compare these results with popula- 
tion size (density and biomass) data and 
the predictions of the thermal equilibrium 
hypothesis. Three sites were chosen on the 
same river, all with similar elevations but dif- 
ferent temperature regimes. Site 1 was lo- 
cated near the upper altitudinal limit for each 
species, whereas sites 2 and 3 were located in 
a warmer, more constant thermal regime 
downstream from a deep-release reservoir. 
Specifically, we sought to test the following 
hypotheses: (1) population size and body size 
will be smaller for each species at site 1 com- 
pared with sites 2 and 3; (2) population size 
will correspond with body size at the same 
sites (i.e., site-specifie ranks of body size and 
population size should be the same); (3) sea- 
sonal patterns of growth will parallel seasonal 
temperature patterns. Body  size/fecundity 
was compared to examine its ability to explain 
among-site differenees in population size. 

Several studies have demonstrated a posi- 
tive correlation between body size and fecun- 
dity in aquatic insects (see Clifford and Boer- 
ger 1974, Kondratieff and Voshell 1950, 
Sweeney and Vannote 1981). Therefore, we 
assume that larger mayflies produee more 
eggs compared with smaller mayflies of the 
same species. We did not attempt to deter- 
mine, as predicted by the TE hypothesis and 
numerous other authors, whether an inerease 
in population size is positively correlated with 
an inerease in population stability. Consider- 
able evidence, however, 
populations are more stable than smaller pop- 
ulations (e.g., smaller populations are more 
susceptible to extinction). 


STUDY SITES AND INSECT 
ALTITUDINAL DISTRIBUTIONS 


The first hypothesis required us to a priori 
rank the study sites according to how closely 
they approximated the optimal temperature 
conditions for the maytlies under investiga- 
tion. Our ranks were based on the altitudinal 
distributions of the inseets. Because the study 
sites were located near their upper limits, the 
unaltered temperature regime at site I was 
considered cooler than the optimum neces- 
sary to maximize body size/fecundity. There- 
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fore, the warmer, dam-impacted sites 2 and 3 
were assumed to be nearer the insects’ opti- 
mal temperature regime. There was no a pri- 
ori reason to separate sites 2 and 3 with re- 
spect to their influence on growth, body size, 
and population size even though they had 
very different temperature regimes. 

The study was conducted in the Upper Col- 
orado River on the western slope of the Rocky 
Mountains in the vicinity of Granby Reser- 
voir, a large (666 km’), deep- release storage 
apodndhene Granby Reservoir is loeated 
38 km northwest of Denver, Colorado. Site | 
was located in a third-order, free-flowing sec- 
tion of the river 4.0 km above the reservoir; 
sites 2 and 3 were loeated 0.4 and 4.0 km, 
respectively, downstream from the dam. Al- 
though differentially influeneed by stream 
regulation, all three sites had similar gradi- 
ents (0.006—0.009), canopy cover, geology, 
riparian vegetation, and elevation (2,593 m, 
2,454 m, and 2,426 m, respectively. (For ecom- 
plete site deseriptions see Rader and Ward 
[1988].) 

Three mayfly species were analyzed in this 
study: Drunella grandis (Eaton), Ephemerella 
infrequens ‘MeDunnough), and Baetis tricau- 
datus (Dodds). Ward (1980, 1986) determined 
the altitudinal distribution of macroinverte- 
brates, including the mavflies of this study, in 
the St. Vrain River, a free-flowing stream run- 
ning from the alpine tundra to the plains on 
the eastern slope of the Rocky Mountains. 
Based on his results, site | of this study (upper 
montane zone) was above the altitudinal dis- 
tribution for D. grandis and very near the 
upper limits for E. infrequens and B. tricau- 
datus. All three species exhibited maximum 
cdlensities at lower elevations in the foothills or 
plains zones. Even though their altitudinal 
upper limits appear to be somewhat higher in 
the Colorado River, probably beeause of its 
larger size compared with the St. Vrain River, 
we concluded that all three study sites were 
located near the upper altitudinal limit for 
each of the three mayfly species. 


METHODS 
Temperature, Growth, and Body Size 


Water temperature was measured continu- 
ously at each site for 1S months using Ryan 
90-day thermographs. Eaeh thermograph was 
checked against a Weksler hand-held ther- 
mometer on a monthly basis and ealibrated 
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prior to placement and following retrieval. 
Daily mean temperatures were used to calcu- 
late annual mean temperatures, annual coefti- 
cient of variation, annual degree days, num- 
ber of days less than 3 C, number of days equal 
to 0 C, length of spring-summer and summer- 
autumn periods, and rate of spring-summer 
increase and summer-autumn decline. 

Annual growth rate analyses and general 
temperature descriptors (e.g., accumulated 
degree days, mean annual temperature, 
etc.) cannot explain site-specific variation in 
aquatic insect body size and fecundity be- 
cause they average over important seasonal 
information. seaconal growth rate analyses 
combined with seasonal temperature profiles 
can, however, provide insights into the rela- 
tionships between temperature, growth, and 
body size. 

Temperature profiles for each site were 
separated into three periods: (1) spring-sum- 
mer, (2) summer-autumn, and (3) winter 
(Fig. 1). The winter period was defined by 
mean daily temperatures <3 C in order to 
inelude the winter warm temperatures at 
site 2. The end of the spring-summer period/ 
beginning of the summer-autumn period was 
set at 15 August, based on temperature peaks 
apparent at sites 1 and 3 (Fig. 1). Therefore, 
the spring-summer period began when the 
mean daily temperature exceeded 3 C and 
ended 15 August. The summer-autumn pe- 
riod began 15 August and ended when the 
mean daily temperature dropped below 3 C. 

Growth, defined as the monthly increase in 
mean biomass of individuals collected per 
sampling date, was determined for D. grandis 
and E. infrequens. Growth for B. Aan nics 

was not analyzed because of difficulty in as- 
signing intermediate-sized instars to the cor- 
rect generation. 

Site-specific differences in seasonal growth 
rates were determined by regressing monthly 
mean mdividual biomass estimates against the 
number of Julian days accumulated over the 
three separate growth periods (spring-sum- 
mer, summer-autumn, and winter). A slopes 
comparison test (analysis of covariance) was 
used to determine among-site differences in 
seasonal growth rates. No transformation was 
necessary because growth was linear over the 
short seasonal periods. 


Population Size and Body Size 


Estimates of population density and _ bio- 
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mass were based on four Surber samples 
(0.09°m each, 240 wm mesh) collected monthly 
cros@ lnenvidth ofthe stueaiwemench sie and 
four artificial substrates. Artificial substrates 
consisted of clay bricks (23 x 19 * 9.5 em) that 
had been in the streanis for one month before 
being sampled. Most of the Surber samples 
enclosed natural substrate particles similar in 
size to the artificial brick substrates. There- 
fore, sampling units from both techniques 
were combined because the sampling areas 
were approximately equal. A simple ¢ test 
indicated that the mayfly population means 
based on Surber samples of natural substrate 
were not significantly (P — .63) different from 
population means based on the artificial sub- 
strate samples. Following identification and 
enumeration, all nymphs were divided into 
0. 1-mm size classes, based on maximum head 
capsule width, and dried at 60 C for 48 hours. 
Mean annual population biomass was deter- 
mined by summing biomass estimates for all 
size classes across all sampling dates. Mean 
annual population density was determined 
after summing the number of individuals in 
each size class at each site. Head-capsule 
measurements were also used to construct 
size-frequency plots for life-history determi- 
nations, including the number of annual 
generations produced. Complete life-history 
information for these species at each site can 
be found in Rader and Ward (1989). 

Female body size (dry weight biomass) of 
late instars was estimated by taking the mean 
of the three largest size classes collected from 
each site. Late mayfly instars have a full eom- 
plement of mature eggs. 


RESULTS 
Temperature, Growth, and Body Size 

A progressive increase in accumulated de- 
gree davs and mean annual temperatures was 
found from site | to site 3 (Table 1). Site 1 was 
characterized by rapid seasonal changes in 
temperature and a long, freezing (0 C) winter 
(Table 1, Fig. 1). In contrast, site 2 demon- 
strated gradual seasonal changes combined 
with winter warm and summer cool tempera- 
tures. Site 3 was intermediate with respect to 
rates of seasonal change but had the greatest 
amount of thermal energy (largest number of 
accumulated degree days, largest maximum 
temperature, and largest annual mean tem- 
perature). 
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Fig. 1. Temperature profiles for each site during [9S1—1982. Individual points represent daily means. (See text for 


explanation of criteria used to determine seasonal separations. Note that the x-axes have been adjusted to facilitate a 
comparison of the length ofeach period at each site. ) 


1990] 


Tas.eE lL. Temperature characteristics of each site. 





Characteristic Site | Site 2 Site 3 





Annual degree days (C) M30 1729" 2052 
Mean annual temperature (C) Soret 8 —O. 
ov. (&) omnes Gd 
Minimum (C) 0.0 18 0.0 
Maximum (C) 18.0 9.8 18.2 
Davs <3 C It) i53° 3101 
Days ee O C loo ed 17 
Rate of spring-summer fT (C/day) 0.14 0.05 0.06 
Length of spring-summer (days) 99 OF 177 
Rate of summer-autuinn | @C/day) 0.20 0.05 0.10 
Length of summer-autuinn (days) 70 115 87% 


Seasonal growth rates were analyzed to ex- 
plain among-site differences in body size for 
each species. The seasonal pattern of growth 
for D. grandis and E. infrequens at site 1 was 
characterized by slow growth during summer- 
autumn, no growth during winter, and a rapid 
increase during the spring-summer period 
(Figs. 2, 3). In contrast, growth at site 2 was 
comparatively fast and continuous through 
each seasonal period including winter. At 
site 3, D. grandis demonstrated a continuous 
growth pattern similar to that observed at 
site 2, whereas E. infrequens demonstrated a 
spring-summer pulsed pattern of growth 
more similar to that at site 1. 

The seasonal growth rate (0.0108 mg/day) 
for the early instars of E. infrequens during 
the summer-autumn period was significantly 
(P = .0001) faster at site 2 than at sites 1 and 3 
(0.004 and 0.003 mg/day, respectively), which 
were not significantly different (Fig. 2). This 
trend continued into the winter period when 
the growth rate at site 2 (0.022 mg/day) was 
again significantly (P = .0001) faster than at 
sites land 3, where growth rate was not difler- 
ent from zero. By the end of winter, the site 2 
population had completed over 76% of its 
growth, and larvae were over five times larger 
than those at sites 1 and 3 (Fig. 2). Therefore, 
the larger body sizes at site 2 can be attributed 
to rapid growth starting at egg hatch and con- 
tinuing through the winter. During the 
spring-summer growth period, individuals at 
sites | and 3 grew significantly (P — .0001) 
faster (0.034 and 0.033 mg/day, respectively) 
than individuals at site 2 (0.014 mg/day). 
However, the body sizes of late instar larvae at 
sites 1 and 3 were still considerably smaller 
than those at site 2. 

The growth rate (0.064 mg/day) of early 
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instars of D. grandis during the siummer- 
autumn period was significantly different 
among sites (P — .OOOL), being greatest at site 
2 and slowest at site | (0.035 mg/day). Winter 
growth was also faster at site 2, with an aver- 
age rate of 0.064 g/day, followed by site 3 
(0.044 mg/day) and then site 1 (0.015 mg/day). 
Winter growth at site | was not significantly 
different from zero. Spring-summer growth at 
site 1, however, was significantly (P ~ .0001) 
faster than at sites 2 or 3, with the fastest 
seasonal rate of increase for this study (0.119 


mg/day). Spring-summer growth rates at 
sites 2 and 3 (0.076 and 0.055 mg/day, respec- 


tively) were not significantly different. 

Bactis tricaudatus was univoltine at site 1, 
but bivoltine with slow and fast seasonal gen- 
erations at sites 2 and 3. In contrast, D. gran- 
dis and E. infrequens had univoltine, slow 
seasonal life cycles at each site. Complete life- 
history data tor each species can be found in 
Rader and Ward (1989). 

As predicted by the first hypothesis, mean 
annual density and biomass of each species 
were much greater at the warmer sites below 
the dam than at site 1 (Table 2). Population 
density and biomass of Baetis tricaudatus and 
E. infrequens were largest at site 2, followed 
by sites 3 and 1. Maximum density and 
biomass of D. grandis were greatest at site 3, 
followed by sites 2 and 1. 

Contrary to the predictions of the second 
hypothesis, dry weights of the largest instars, 
and population sizes of each species, did not 
correspond when ranked by sites (Table 2). 
Although the largest late instars of E. infre- 
quens occurred at site 2 2. where the density 
was greatest, the largest instars of B. tricauda- 
tus aad D. grandis did not oceur in the largest 
population. Body size and population size cor- 
responded in only two other instances; the 
smallest late instars of B. tricaudatus and the 
intermediate-sized late instars of D. grandis 
occurred in the small- and arenmedintes sized 
populations, respectively, at site 1. 

The largest B. tricaudatus larvae occurred 
at site 3, whereas its maximum population 
density and biomass occurred at site 2. Late 
B. tricaudatus instars at sites 2 and 3 were 
over two times larger than late instars at site 1 
(Table 2). The largest-sized larvae of D. gran- 
dis oceurred at site 2. Site 3, which had the 
largest population of D. grandis, had_ the 
ret late instars (Table 2). 
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Fig. 2. Seasonal growth rate patterns for Ephemerella infrequens, Site Lis represented by a solid line and diamonds, 
sites 2 and 3 by a dotted line and circles, and a dashed line and triangles, respectively. Each symbol (diamond, circle, 
and triangle) represents the mean biomass for a single sample. 
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Fig. 3. Seasonal growth rate patterns for Drunella grandis. (For further description, see Fig. 2 legend.) 
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TABLE 2. Mean annual population density (#’sm ~)and biomass (mgm ~), plus mean individual size (mg dry wt.) for 
late instars of each species. Values in parentheses for the population parameters are the percentage of mean 
represented by the standard error. Values in parentheses for size estimates indicate the number of individuals used to 
determine each mean. 











Site ] mite 2 Sitters 
Population size Body size Population size Body size Population size — Body size 
Species Density Biomass Density Biomass Density Biomass 
B. tricaudatus Tot 142.6 0.730 7720 EO Ul rerenl) 4015 1162.8 1.960 
(21%) (17%) (311) (18%) (21%) (33) (1S%) (27%) (35) 
E. infrequens 367 902.6 3.500 1034 4964.7 4.930 644 1188.5 2.800 
(25%) (23%) 23h (23%) (25%) (35) (24%) (26%) (25) 
D. grandis 102 [580.9 14.930 SO 2368.4 18.110 215 3300.4 14.220 
(31%) (18%) (16) (10%) (22%) (9) (14%) (20%) (17) 
DISCUSSION cially applicable to cool-adapted boreal spe- 


cies (see Edmunds 1982). 

Early instars at site 2 began growth in much 

The winter warm and summer cool condi- cooler summer-autumn temperatures; meta- 
tions of site 2 allowed rapid continuous growth — bolic costs were low and growth rates fast 
of E. infrequens and D. grandis, which pro- compared to those at sites 1 and 3. Winter 
duced larger instars and greater fecundity. temperatures, which varied slightly around 
Site-specific explanations of growth patterns 2 C, were not sufficiently cold to inhibit 
and body-size differences for D. grandis and growth, which continued at a rapid pace. 
E. infrequens are consistent with the TE hy- | Growth appeared near completion before the 
pothesis. Vannote and Sweeney (1980) pro- vernal rise in temperature, thus leaving 
posed that the seasonal pattern of growth for — plenty of time for maturation and emergence. 
aquatic insects, as for other small ectotherms The rapid completion of growth probably re- 
(e.g., Phillipson 1981), may be determined by — sulted in the extended emergence of E. infre- 
the interaction between temperature and — quens and the addition ofa second generation 
body size. Smaller instars, which have a large — of B. tricaudatus at site 2. 
surface-to-volume ratio, will have a higher The early instars of D. grandis (July) and E, 
metabolic rate than larger instars at the same — infrequens (August and September) began 
temperature. At site 1 the smallest D. grandis — growth at site 3 during the warmest coe of 
and E. infrequens instars appeared during the — the year. August and September were, on 
warmest months of the year (July, August, average, 7-S C warmer at site 3 than site 2. 
and September). High summer-autumn tem- Although E. infrequens did not grow, D. 
peratures coincident with small instars at grandis early instars grew rapidly during the 
site 1 likely resulted in large metabolic costs | summer-autumn growth period. Because D. 
and, therefore, slow growth rates and possibly — grandis early instars were approximately two 
high mortality rates. Both species ceased to — times larger than E. infrequens early instars, 
grow during the freezing winter tempera- they probably had lower metabolic costs. This 
tures. In the spring, winter survivors experi- allowed them to grow at the warm summer- 
enced a rapid increase in temperature and autumn temperatures. Although both species 
thus a relatively short period (99 days or less) stopped growing during winter at site 1 but 
to complete growth and maturation. All else — continued to grow during winter at site 2, only 
being equal, the magnitude and length of | D. grandis continued winter growth at site 3. 
summer-autuinn temperatures when coinci- The fact that D. grandis was over three times 
dent with early instars, plus the rate of vernal — larger than E. infrequens at the beginning of 
rise, limit growth and body size/fecundity and — winter may explain its ability to grow in winter 
probably have an important influence on the — at site 3 in contrast to E. infrequens. 
geographic distribution and upper altitudinal These site-specific explanations of growth 
limits of aquatic insects. This may be espe- are consistent with the TE hypothesis 


Temperature, Growth, and Body Size 
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suggesting that growth rate and consequently 
body size and fecundity are determined by 
the length of time individuals are exposed to a 
specific optimal range of temperatures. Other 
factors, however, that may also influence 
maylly growth rates (e.g., food abundance; 
Sweeney et al. 1986) were altered by the ef- 
feets of stream regulation. For example, con- 


stant flow conditions and the addition of 


planktonic diatoms from the reservoir en- 
hanced food quality and quantity at sites 2 and 
3 (Rader and Ward 1989). 

If summer temperatures increase meta- 
bolic costs, causing growth to slow or stop, 
then the rate of vernal rise and autumn de- 
cline determines the amount of time individu- 
als are exposed to optimal temperatures and, 
therefore, the amount of time available for 
growth. Growth of both D. grandis and E. 
infrequens continued as long as temperatures 
remained between 2 C and 10 C. However, 
when temperatures exceeded this range, 
growth slowed or stopped. Where growth was 
continuous (site 2), temperatures were always 
within this range. The optimal temperature 
range for these two species appears to lie be- 
tween 2 C and 10C. 


Population Size and Body Size 


Winter warm and summer cool conditions 
at site 2 and the long spring-summer period at 
site 3 probably accounted for the multivoltine 
life cycle of B. tricaudatus at these sites. Stan- 
ley and Short (1988) suggested that population 
size may remain unaltered or even increase in 
warmer than optimal conditions if faster 
growth rates and shorter generation times 
compensate for smaller body sizes and lower 
fecundity. For some aquatic insects, warmer 
than optimal temperatures may offer a trade- 
off between body size/fecundity and genera- 
tion time. Will they maximize reproductive 
effort by producing fewer, larger individuals 
(slow growth and a univoltine life cycle) or 
many, smaller individuals (fast growth and 
multivoltine life cycle)? These data demon- 
strated that Baetis may have the genetic plas- 
ticity necessary to respond to such tradeoffs. 
As temperatures approached optimality 
(site 2) from cooler conditions (site 1), both 
voltinism and body size increased within the 
same population. When comparing site | with 
the much warmer and very different tempera- 
ture regimes downstream from the reservoir, 
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we found that our data support the predictions 
of the TE hypothesis. The largest bedy size 
and population size of each species occurred 
downstream from the reservoir. However, 
site-specific comparisons of body size and 
population size did not correspond as pre- 
dicted by the TE hypothesis. The therimal 
regimes at sites 2 and 3 were apparently snffi- 
ciently similar that temperature did not have 
an overriding influence on population sizes. 
The TE hypothesis assumes that higher fecun- 
dity is equivalent to larger population size. 
Sources of mortality at every stage of develop- 
ment, eggs, nymphs, and adults (Sweeney 
and Vannote 1982, Butler 1984, Peckarsky 
1984, Gilliam et al. 1989), which may vary 
across sites, may interrupt the translation of 
higher fecundity into higher population den- 
sity and biomass. Numerous factors, in addi- 
tion to the influence of temperature on body 
size/fecundity, will undoubtedly influence 
the geographic or altitudinal variation in 
aquatic insect population size and stability. 
Toward the center, and probably over most of 
a species range, other sources of mortality and 
determinants of reproductive success should 
have a greater influence on aquatic insect pop- 
ulation size. Temperature at the edge of a 
species range may, however, be more limit- 
ing, compared to other factors, in determin- 
ing the extent of a species distribution. 
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